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We present an analytical treatment of coherent excitation of a Two-Level Atom driven by a 
far-off resonant classical field. A class of pulse envelope is obtained for which this problem is 
exactly solvable. The solutions are given in terms of Heun function which is a generalization of 
the Hypergeometric function. The degeneracy of Heun to Hypergeometric equation can give all 
the exactly solvable pulse shapes of Gauss Hypergeometric form, from the generalized pulse shape 
obtained here. We discuss the application of the results obtained to the generation of XUV. 
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I. INTRODUCTION 

The two-level system (TLS)[l-3] is a very rich and 
useful model that helps to understand physics of many 
problems ranging from interaction with electromagnetic 
fields to level-crossing[4-6]. For example, interaction 
of a beam of atoms in Stcrn-Gcrlach apparatus [7], and 
Bloch-Siegert shift [8] can be understood using TLS. Re- 
cently TLS has been extensively studied as a quantum 
bit (qubit) for quantum information theory [9]. Two-level 
atom (TLA) description is valid if the two atomic levels 
involved are resonant or nearly resonant with the driv- 
ing field, while all other levels are highly detuned. TLS 
can be realized exactly for a spin-1/2 system, and, ap- 
proximately, for a multi-level system in a magnetic field 
when all other magnetic sub-levels are detuned far-off 
resonance. 

When the frequency of the driving field is in resonance 
with the atomic transition frequency, the Schrodingcr 
equation for the time evolution of state amplitudes is ex- 
actly solvable for any time dependence of the field Q(t). 
The transition probability is given as 

/oo 
n(t)dt. (i) 

Here A is the area of the pulse envelope. Interestingly 
this transition probability vanishes when A is an even 
integer multiple of n (CPR). For odd integer multiple 
of 7r we get complete population inversion (CPI) while 
half-integer multiple of 7r gives equal coherent superpo- 
sition of the initial and the final states. Several exactly 
solvable models for the TLS have been proposed in the 
past [10-24] where solutions to the Schrddinger equation 
are expressed in terms of known functions like Hyper- 
geomteric functions. Several approximate solutions have 
also been proposed based on perturbation theory and the 
adiabatic approximation[25,26]. 

Recently, the topic has been in a focus of research re- 
lated to generation of short wavelength radiation[27,28]. 
A two-level atomic system under the action of a far-off 
resonance strong pulse of laser radiation has been con- 



sidered and it has been shown that such pulses can ex- 
cite remarkable coherence on high frequency far-detuned 
transitions; and this coherence can be used for efficient 
generation of UV and soft X-ray (XUV) radiation [28]. 

To describe excited coherence, we are interested to un- 
derstand the mechanism of breaking adiabaticity that 
leads to excited coherence in the system when the laser 
pulse has already passed. Thus we are interested going 
beyond classical electrodynamics [29]. Indeed, an electric 
field causes polarization of dielectrics is given by 

P(t,r)= f dt' X (t-t')£(t',r), (2) 

J — oo 

where x(r) is the dielectric response function. It is im- 
portant to note that once the field is removed, the polar- 
ization adiabatically returns to practically zero. Break- 
ing of adiabaticity is especially difficult when the fre- 
quency of the applied field is far from the atomic res- 
onance. Finding exact analytical solutions for such a 
problem will not only supplement numerical simulations 
but will also be useful in understanding the underlying 
physics. 

In this paper, using a proper variable transformation, 
we find a class of pulse Sl(t) for which the Schrodingcr 
equation for the time evolution of the state amplitudes 
can be transformed into the well known Heun equation. 
The solutions are given in terms of the Heun function 
which is a generalization of the Hypergeometric function. 
Using the degeneracy of Heun to Hypergeometric equa- 
tion, Bambini-Bcrman model[21] can be generalized to 
this model. 

The paper is organized as follows. In section 2, we 
briefly describe our system and obtain the equation of 
motion for the state amplitudes. In section 3, we present 
the exact solution of the problem in terms of the local 
Heun solutions HL It is well established that the Heun 
equation reduces to the Gauss Hypergeometric equation 
in several ways so we discuss this degeneracy briefly. We 
also discuss one of the confluent cases of the Heun Equa- 
tion i.e the Confluent Heun Equation and find the exact 
solutions. In section 4, we give some specific examples 
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of the pulses for which we have found solutions. Pulse 
shapes are asymmetric in time except the Rosen-Zener 
pulse. The Hyperbolic secant (Rosen-Zener Model), gen- 
eralized Rosen-Zener (Bambini-Berman) Model are in- 
cluded in this class f2(£) as a special case. We also give a 
new model for a Smooth Box Pulse which takes care of 
non-analyticity at the edges by introducing a parameter 
5. By modulating this parameter we can modulate the 
box width. In section 5, we discuss the application of the 
results obtained here to the generation of XUV. We also 
estimate the level of the XUV field that can be generated 
by using the excited coherence. 



II. TWO-LEVEL ATOM: EQUATION OF 
MOTION 

The equation of motion for the probability amplitudes 
for the states \a) and \b) (see Figl(a)) of a Two Level 
Atom (TLA) interacting with a classical field is given as 



C h = 



h 



cos{vt)e- iut C a , 



(3a) 
(3b) 



where hu> is the energy difference between two levels, 
p is the atomic dipole moment; E(t) = £(t)cosvt (see 
Figl(b)). In the Rotating Wave Approximation (RWA) 
we replace cos(z4)e ±Jwt — > e ±iA /2 where A = lu — v, is de- 
tuning from resonance. Introducing £7(£) = p£ (i)/27i[30], 
Eq.(|3]) reduces to 



C a = ifl(t)e iAt C b 



C b = in*(t)e- iAt C a 



(4a) 
(4b) 



which have an integral of motion \C a \ 2 + \C b \ 2 = 1 [31] . 
There are a variety of ways to approach the problem 
of solving for C a (t) . One method is to define f(t) = 
C a (t)/C b {t). For the function /(*), Eq.Q yields the fol- 
lowing Riccati Equation [28] 



/ + ifT (t)e- lA 7 2 - iSl(t)e tm = 



(5) 
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FIG. 1. (Color online) (a)Two-level atomic system, atomic 
transition frequency uj — ui a , — uj b , detuning A = ui — v and 
Rabi frequency Q(t) — p£(t)/2h. (b) Classical electromag- 
netic field E(t)= sech(at)cos(i4) 



Then \C a (t\ = \f(t)\/y/l + \f(t)\ 2 . Alternatively, we can 
get a second order linear differential equation for C a (i), 
from Eq.Q 



C a (t) 



c a (t) + \n\ 2 c a (t) = o. 



(6) 



The general solution for Eq(|6]) has not been found yet, 
however there are solutions for several cases in terms of 
special functions. To find a solution for Eq.Q we intro- 
duce a new variable 



(7) 



subject to the condition that tp(r) is real, positive and 
monotonic function of r and ipo < (f < ip\. In terms of 
the variable ip and the dimensionless parameters 



A 

a 



a 



at, f) = — , 7 
one may write Eq.(|6]), for real £(r) in the form 



(8) 



C„ 



2 £-2 



c„ 



C a = 0, (9) 



where a prime indicates differentiation with respect to ip 
and O(t) = 7^(t). Let us determine the condition under 
which Eq.([9]) has the form 



c>) + P{<p)dM + Q(<p)c a (<p) = o. 

Using Eq.( |9|10p and some trivial alebra we get, 




(10) 



(11) 



III. HEUN EQUATION 

Bambini-Berman studied the case in which Eq.( 10]) has 
the form of a Gauss Hypergeometric equation which in- 
cludes Rosen-Zener Model as a special case. Now let us 
consider when Eq.(10) is of the form of Hcun equation 



[32,33] with the independent variable ip. 



d 2 Cg 

dip 2 



dC a (abip - q)C a 



tp ip — 1 ip — c) dip ip{ip — l)((p — c) 



= 0, 



(12) 

where a,b,c,q,u,v,w are parameters with c ^ 0, l.(c > 1). 
The parameters are constrained, by the general theory of 
Fuchsian equations, as 



u + v + w = a + b+l. 



(13) 



From Eq.(12| and Eq.(|10| and some algebra we get 
(p = 2<p(l - ip)/((j, + \(p). 



(14) 
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Equivalently the parameters of the Heun Equation 
Eq.(12| are given as 



i/3/J. 
2 ' 



1 i/3(X + fi) 



2 

a = 0, 



1 

2' 



, 1 0X 



(15) 



For tp(r) to be a monotonically increasing function of r, tp 
must be real and positive i.e /i > 0, A//i > — 1. The time 
variable r as a function of tp is obtained by integrating 
Eq.(14| which gives, 



2T = ln[^/(l-(^)^ +A ] 



(16) 



The general solution for Eq.(12), which has regular sin- 



gularity at tp = is given in terms of the Heun local 
solutions, HI ((f) as, 

C a = P 1 ^ 1 -"ffl[c, q + (1 - u)((c - l)v + a + b - u + 1); 
a — u + 1, b — u + 1, 2 — u, v; tp] + P2HZ [c, g; a, b, u, v; tp] , 

(17) 

where the constants, Pi, P2 can be found using the initial 
conditions of the system. In the limit r — > 00, the popu- 
lation left in the level \a) can be obtained by substituting 
p — > 1 in Eq. ( p~7| ) . The form of the pulse can be obtained 
by equating Eq.Q and Eq.(12) which gives 



n(r) = 



Atp(l - p)(abp - q) 
(c-tp) 



1/2 



1 



pL + Xp 



(18) 



where tp(r) is given by Eq.(16)[34]. In Fig (2) we have 
plotted the pulse envelopes' of the classical field, given 
by Eq.( 18 ), for which the two-level atom problem can be 



exactly solved. They also show the effect of the asym- 
metric parameters A and ab respectively, for p, = 1, on 
the symmetry of the shapes. Pulse shapes showing the 
effects of other parameters can also be plotted easily from 
Eq([l8]). 



There are three kinds of solutions to the Heun equa- 
tion Eq.(fl2l). Local Solutions HI, Heun functions H/ and 



Heun Polynomials Hp[35-37]. The series solution Eq.( 17) 
is written as [33] 



HZ [c, q; a, b, u, v; p] = s : 
3=0 



J — 



1 + — (p- 
uc 



J=2 



(19) 



where Sj obeys the three term recursion relation 

(i - 1 + a)(j - 1 + 6)sj_i - - 1 + u)(l + c) + vc 
+ a + b+ l- u — v]+ q}sj + (j + l)(j + u)s j+ i = 0, 

(20) 



with the initial conditions 



S = 1, Si 



q_ 

uc 



and Sj = 0, if j < 0. (21) 



The solution Eq.( 19 1 is valid only within a circle centered 
at the origin tp = whose radius is the distance from the 
origin to the nearest singularity ip = 1 or tp — c. For 



c > 1, the radius of convergence is 1[33]. From Eq.(20l, 
we can say that Heun function remains the same with 
the exchange of the parameters a and b. 



III.l. Degeneracy to the Hypergeometric Models 



It can be easily verified that the Heun equation Eq.( 12 ) 



can be reduced to the Hypergeometric equation in several 
ways [33]. They are 



(22a) 
(22b) 
(22c) 

Let us now consider the simplest case of c = 0, q = 0. 
Then for a + b = and 1/2-v = -0/2, Eq.Q reduces 



C = 1, q = ab, 
w = 0, q = cab, 
c = 0, q = 0. 



to standard form of the Gauss Hypergeometric equation 



d 2 Cg 

dp 2 



-(1 



b)p 



tp(l - ip) 



dCg 

dp 



abCg 
tp(l - tp) 



0. (23) 



where r = 1/2 — 0/2. The general solution for Eq.(23) 
is 



C a = Pi^- r F[6 _ r + l,o-r + l;2- 
+ P 2 F[a,b;r;p}, 



r;tp] 



(24) 



where the constants, Pi, P2 can be found using the initial 
conditions of the problem. We write the hypergeometric 
series F( 2 ,i) [a, b; c; p] as F[o, b; c; tp]. The population left 
in the state \a) is given as 

C af = PiF[6-r+l, a-r+l; 2-r; l]+P 2 F[a, b; r; 1]. (25) 

Subsequently if (a + b) = Xi(3 and v — 1/2 — (a + b) = fiif}, 
we have the generalized Rosen-Zener Model as discussed 
by Bambini and Berman[21]. One can summarize the de- 
generacy of the Heun to Hypergeometric model as follows 



HI [1, ab; a, b, u, v; p] — F[a, b; u; tp], (26a) 
HI [c, cab; a, b, u, a + b — u + 1; tp] = F[a, b; u; tp], (26b) 
HI [0, 0; a, b, u, v; p] = F[a, b; a + b - v + l;tp]. (26c) 



III. 2. Confluent Heun Equation 

The Confluent Heun Equation is one of the four conflu- 
ent forms of Heun's equation which is obtained by merg- 
ing the singularity at tp = c that at tp = 00. Now we 
have a regular singularity at p — 0,1 and an irregular 
singularity at tp = 00. In this paper we will consider the 
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FIG. 2. (Color online) Pulse shapes given by Eq(18). (a) Pulse shapes with varying A and c = 2, q = — l,ab — 0. (b) Pulse 
shapes with varying c and A = 2, q = — 1, ab = 0. 




FIG. 3. (Color online) Pulse shapes given by Eq(31).(a) Pulse shapes with varying A and p = — q = 1. (b) Pulse shapes with 
varying q and A = 2,p = 0. 



following non-symmetrical form of the Confluent Heun 
equation: 



d 2 C a 



dC a 



pip- 



-C a = 0. (27) 



dip 2 ' \ip ' ip — I J dp ' ip(ip — 1) 
Similar to the Heun case, we have the same differential 



equation for ip i.e Eq(14|. For the Confluent Heun Equa- 
tion, the possible values of the asymmetric parameters 
are 



i/3/j, 

~2~' 
i/3/j, 

~2~' 



V = n > P 



-3) 



(28a) 



1 

7J = - 

2 



ij0(A + m) 



p = (28b) 



The general solution of the Confluent Heun Equation 
Eq.(27l is given as 



C a = PiHZ (c) [0, u - 1, v - l,p, q + (1 - uv)/2, p] 

+ V 2 p 1 ~ u m^ [0, 1 - u, v - l,p t q + (1 - uv)/2, ip], 

(29) 



where Pi, P2 can be found using the initial condition of 
the system. It is worth mentioning here that, the general 
solution to the Gauss Hypergeometric differential equa- 
tion Eq.( 23 ) can be expressed in terms of the Heun func- 
tions Hl^ 1 as 

C a = Pi(</3 - 1)- Q ffl (c) [0, a - b, -1 + r, 0, ((r - 2a)6 - r 
+ ra + l)/2, 1/(1 - ip)} + P 2 (ip - l)- b Rl^[0, b - a, 
- 1 + r, 0, ((r - 2a)b - r + ra + l)/2, 1/(1 - p)]. 

(30) 

The form of the pulse can be obtained by equating Eq. ^ 
and Eq.(27) which gives, 



f%) 



[Ap{p-l){pp + q)] 1/2 
fi + Xip 



(31) 



where p(r) is given by Eq.(16)[34]. The constraint of A 
and fi is also the same as for the Heun case discussed 
earlier. Fig. (3) shows the pulse shapes for which the 
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two-level atom can be reduced to the Confluent Heun 
equation. It also qualitatively shows the effect of the 
asymmetric parameters p and q on the symmetry of the 
pulse shapes. A = corresponds to the symmetric pulse. 



IV. SOME EXAMPLES 

In this section we will consider some specific exam- 
ples of pulses corresponding to Heun and Confluent Heun 
equations. Interestingly we will also find a better ap- 
proximation for a box pulse by introducing a parameter 
S which takes care of non-analyticity of the pulse at the 
edges. 



IV.l. fi 4 (t) = Q sech(at) / s/8 - tanh(crf), 8 > 1 

For this pulse, using the scaling parameters Eq.Q, 
Eq.© gives 



C a (r)- 







1 ( 1 - 2<Stanhr + tanh r 



2 \_ 5 — tanhr 
7 2 sech 2 r 



Ca{r) 



(32) 



-C a (r) = 0. 



5 — tanhr 
Let us now define a new variable as 

I + tanhr 
VKV = o • 



(33) 



In terms of the variable ip, Eq.(32) reduces to the Heun 
equation 



c a - 

where 



ip ip — 1 ip — c 



a 



abtp — q 



tp((p - l)((f - c) 



it 



v 



1 iB 
— h — 

2 2' 



1 iB 

2 ~ ~2~' 

7 2 1 
g=-y, o = 0, b=- c 



2' 
8 + 1 



C a = 0, 
(34) 

(35a) 
(35b) 




— 8=1.000 
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FIG. 4. (Color online) (a) Pulse shapes for different value of 
8. (b) The time dependence of the population in the state \a) 
for Q.s{t) pulse for different values of 8 > 1. For calculation 
we take a — 0.0Slj c , varying 8. 



Now for this pulse, using the scaling parameters Eq.p]), 
Eq.© gives 



c a (ry 



+ - (1 - tanhr) 



In terms of the variable ip, Eq.(39) reduces to 

v 



where, 



p p-1 



a 



C* a (r)+ 7 2 (l+tanhr)C a (r) - 0. 

(39) 

(40) 



C a = 0, 



1 iB iB 
u = , v = 1 H 

2 2' 2 



The general solution to Eq. ( 40 1 is 



C{p) a = Pi(p - 1)«F[C, $ - 1 + u + w; u; 49] 
P 2 p 1 - U (p - 1)*F[£ + «, £ + 1 - «; 2 - u; p], 



(41) 



(42) 



where, 



I - v 



1-v 



(43) 



From Eq.([32D we see as r -» -00, y> -> and r -» 00, an( j are given be Eq.flilj). Using the initial condi 



ip — > 1. The initial conditions for our system are 

C a (r^-oo) = 0, |C 6 (r -oo)| = 1. (36) 



tions Eq. ( 36 ) we get Pi — and 

o _ 7 



The complete solution to Eq.(34), satisfying the initial 
conditions Eq([36|), is 



V2(u- 1)(-1)K+V3) 



(44) 



C a (p) = ^ y -»W\c, q + (1 - u)((c - 1)« 
+ a + 6 — u + 1); 6 — tt + 1, — u + 1,2 — u, v,(p]. 



where a, 6, c, g, u, u, ti) are given be Eq.(35 ). Let now con- 
sider a case in which 6=1. So the pu 



f2i(t) = ilov 7 ! +tanhcrf 



se has the form 
(38) 



Fig. (4) shows the plot of population in the state \a) corre- 
sponding to the pulse Q$ satisfying the initial condition. 



IV.2. fi+(i) = Q sechat(y/1 + tanhat) 

For this pulse, using the scaling parameters Eq.©, 
Eq.© gives 



6 



Ca(r) 



i/3 + -(1 - 3tanhr) 



C a (r) 



(45) 



+ 7 2 sech 2 T(l + tanhr)C Q (r) = 0. 



In terms of the new variable ip, Eq.(45l reduces to the 
Confluent Heun equation. 



p-l 



C a =0, 



(46) 




where, 



W 1 iP 9 2 

u = -y, «= 2 + y ; a = -27. 



(47) 



The complete solution to Eq.(46) satisfying the initial 
conditions Eq.prjl) is 



C a (ip) = (|rqU 1 + z/3/2,-1/2 + ifi/2, 



2 1 2 ,l/2-f3 2 /8-tf3/8,<p]. 



(48) 



FIG. 5. (Color online) (a) Pulse shapes for Q±(t) = 
r2osechat(\/l ± tanhat). (b) Time dependence of population 
in the state |o) for the Pulse shapes in(a). In calculation we 
take flo = 0.02uj c , a = 0.08o; c , A = 0.2u; c 



IV.4. Smooth Box Pulse 

One of the simplest and exactly solvable pulse shapes 
is a Box Pulse. Indeed it is a non-analytical pulse but 
it gives information about the basic oscillatory nature of 
solution (probability amplitude). Let us define our pulse 
as 



IV. 3. Q_(t) = fi sechcrf(vT - tanhai) 

For this pulse, using the scaling transformation Eq.Q, 
Eq.© gives 



C a (r) 



i/3 (1 + 3tanhr) 



Ca(r) 



(49) 



+ 7 2 sech 2 r(l - tanhr)C a (r) = 0. 



In terms of the new variable <p, Eq.(49l reduces to the 
Confluent Heun equation. 



C„ 



where, 



tp p-l 



1 

2 ~ ~2' 



c' a + ^C a = 0, 



V = 2 7 • 



(50) 



conditions Eq.(|36[), is 



n(t) = n e(t)Q(t Q - 1), t >o 



(53) 



where, Q(t) is a unit step function. The solution for 
Eq.© corresponding to the box pulse is 



C a (t) 



v/AVT 



5 i(A/2)t sin ^ A 2/ 4 + 2) i; t < tQ 

(54) 

The oscillatory nature of the solution \C(t)\ is evident 
from the sine function. Let us consider the pulse shape 
of the form 



n,(t) 



f^osechat 
\/8 — tanhcrf 



S = 2c - 1 



(55) 



where c is one of the singularities of the Heun Equation. 
Assuming c > 1 gives 8 > 1 . A pulse shape of the form 
Eq.(55) is positive definite and it vanishes at r = ±oo. 



l/2-2j 2 -(3 2 /8-if3/8,p] 



(51) 
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-5 



(52) 
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In Fig(5) we have plotted the pulse shapes fi±(r) and 
the corresponding time evolution of the probability am- 
plitude for state \a). 



FIG. 6. (a) Box Pulse for 5-1 = 10" 9 . (b) Time dependence 
of population in the state \a) for the Box Pulse Qs(t). In 
calculation we take fio = 0.02oj c , a — 0.08cj c , A = 0.2lu c 



7 



Let us see what happens when 5 approaches but never 
reaches to 1. We see from Fig(4a), that as S approaches 
to 1, the pulse become more and more broad there by 
making it a better approximation for a box pulse (taking 
care of non-analyticity at the edges). The general solu- 
tion for the pulse of the form Eq.(55 ), is given by Eq.( 17) 



where the asymmetric parameters are given by Eq.(35) 



V. DISCUSSION 

The obtained results can be applied to the generation 
of X-ray and UV (XUV) radiation, which is one of the 
main topics in modern optoelectronics and photonics |38j . 
Recent progress in ultrashort, e.g. attosecond, laser tech- 
nology allows searchers to obtain ultra-strong fields [55] , 
Interaction of such strong and broadband fields with a 
two-level atomic system, even under the action of a far-off 
resonance laser radiation is of current interest [2"TII2"51,|4"0T - 
132] . Strong short laser pulses can excite remarkable co- 
herence on high frequency transitions; and this coherence 
can be used for surprisingly efficient generation of XUV 
radiation [27] HH1 SI1 S3 ■ In the first step we excite the 
atoms (e.g., from the Is to 2s states of or He + , etc.) via 
a short pulse of femto- or attosecond radiation e.g., from 
a conventional Ti-sapphire laser system). The excitation 
occurs due to the coherent coupling between Is and 2p 
and then 2p and 2s. In the second step, we apply another 
pulse which scatters off the Raman coherence (prepared 
in the first step) , generating short wavelength anti-Stoke 
radiation as depicted in Fig(7). The generation of radia- 
tion is a coherent process that (contrary to conventional 
superfluorescence) does not require population inversion 
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FIG. 7. (Color online) Field configuration and level structure 
of H or He + . All population is initially in the ground state \b). 
First, the strong short far-off resonant pulse with frequency 
v\ is applied to the system to excite coherence between levels 
\b) = | Is) and |c) = 1 2s) , and then the second pulse with 
the frequency ^3, which is close to the transition between 
levels 2s and 2p, is applied to generate XUV pulse with higher 
frequency va 



-•— Ic> 



-|a>\ ; 



-Id) 
Q 




(a) 



V (b) 



FIG. 8. (Color online) Two-stage generation scheme for X-ray 
generation, (a) Applying a strong pulse allows one to excite 
an atomic system by transferring population to electronic ex- 
cited states, (b) Coherence is then induced by applying a 
resonant field. 



(see Appendix) . The higher efficiency of coherent process 
has been demonstrated in various spectral regions [43-50]. 

We have analytically calculated above that the level 
of excited coherence when a two level atom is driven by 
a ultra-short intense pulse. The coherence is sufficiently 
large that this can be used for nonlinear generation of 
XUV radiation, i.e, see Figs(4b, 5b, 6b), coherence can 
be of the order of 0.1. It is instructive to estimate the 
level of XUV field that can be generated by using this co- 
herence. After an ultra-strong and short pulse, we apply 
a strong resonant and relatively long pulse. The applied 
probe pulse ^3 and generated signal are coupled to 
each other via coherence excited in the medium (Rabi 
frequencies are defined as ^3,4 = ps^E^^/K). Hence, 
the propagation equation for £l± is given by 



9^4 

dz 



(56) 



where p ab is the appropriate atomic coherence (see 
Fig. (8)), and 774 = k^p\N/2U, where p 4 is the dipole 
moment at the transition between . The corresponding 
equation for the density matrix coherence p ab is 



Pab = -TabPab + i^i(Paa ~ Pbb) ~ «^3Pcb, 

and, for short pulses, p ab ~ — iQ,^Tp cb . Then, 
estimate the intensity of the signal field, by 



n 4 = 



kiLp 2 ab N 
2h 



Pcb^3T, 



(57) 



we can 



(58) 



where k 4 is the wavenumber for signal radiation, L is the 
length of the active medium, p ab is the dipole moment 
at the transition between a and b levels, r is the time 
duration of the pump laser pulse. Using the parameters 



N ~ 

n 3 r 



10 16- 

= 1- 



19 
10 3 



cm , p ab ~ ID, L = 100 pm, p cb = 10 
= 1 ps, A = 10 nm, we obtain energy : 



10 nJ — 1 pJ, This estimate shows the promise of the 
approach. This estimate is valid on the time scale when 



the collisions in the plasma destroy the coherence. It 
occurs at the times of order St — 1/acN ~ 1 ps, where 
a is the atomic cross-section for atomic collisions that 
destroy the excited coherence. 



VI. CONCLUSION 

In this paper we have found several analytical solutions 
for a two-level atomic system under the action of a far-off 
resonance strong pulse of laser radiation. The solutions 
are given in terms of Hcun function which is a gener- 
alization of the Hypergeometric function. The Rosen- 
Zener and Bambini-Bcrman Model belongs to this class 
of pulses as special cases. A better approximation for 
box pulse is also obtained here which take care of non- 
analyticity at the edges by introducing a parameter 5. 
The results obtained here have applications to the gen- 
eration of XUV radiation and the estimate shown above 
shows a good potential for a source of coherent radiation. 
The technique used here to get the exactly solvable pulse 
shapes can be generalized to appropriate time dependent 
detuning A = A(t) cases and produce more exactly solv- 
able models (to be reported elsewhere) [51] 
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Appendix A: Generation of radiation by a two-level 
atomic medium with excited coherence 



Let us assume that a two-level atom has some small 
initial coherence p° ab = a/ p aa Pbb- Note that in this pa- 
per, we consider the case when there is no population 
inversion, p° aa < p bb . The density matrix equations for 
atomic coherence are 



dpab 
dt 



iVt(p aa - Pbb), and 



(Al) 



8_ 

Ft 



(Paa - Pbb) = -2iflp ab . 



(A2) 



Then, for the retarded frame 

r = t--, (A4) 
c 

the propagation equation for a resonant field is given by 

-iVPab, (A5) 



dz 



where rj = 3X 2 N'y / (8tt) is the coupling constant. Intro- 
ducing 



= 2 



n dt, 



(A6) 



Eq. ( A5 ) can be rewritten as 
8 2 9 



= — rj sm( 



dzdr 

where 6 can be determined from initial condition as 



(A7) 



Solution of Eq. ( A7 1 is given by 



9 = cf>[l- Jb(2VP=F)], 
and the Rabi frequency is 

Irjz 



Q = 4>Ji(2y/rjzT)^ 



(A8) 
(A9) 

(A10) 



The energy of the generated short wavelength pulse can 
be calculated as 



c 

4tt 



/oo 
\E\ 2 dt = Az N paa huj a b, (All) 
-oo 



and it is equal to the energy stored in the medium after 
excitation. Also it is important to note that the absence 
of population inversion does not influence much of pulse 
energy because of coherent interaction of the radiation 
field with the atomic medium. 

The time duration of the generated pulse is of the order 

of 



47T 



.pulse- mX 2 z ^ 

and it gives the power of the pulse be 
\ 2 zN „ 

Ppulse = —^—AzN~/ r p aa hUJab, 



(A12) 



(A13) 



the solution (by neglecting relaxation processes) is 



X 2 zN 

where the factor shows the brightness of the source 

_ 47T 

in comparison with spontaneous emission of incoherent 
source. 



Pab = iPab^nl 



(A3) 
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